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Abstract 

Memory dephasing and its impact on the rate of entanglement generation in quantum repeaters 
is addressed. For systems that rely on probabilistic schemes for entanglement distribution and 
connection, we estimate the maximum achievable rate per employed memory for our optimized 
partial nesting protocol, when a large number of memories are being used in each node. The above 
rate scales polynomially with distance, L, if quantum memories with infinitely long coherence times 
are available or if we employ a fully fault-tolerant scheme. For memories with finite coherence 
times and no fault-tolerant protection, the above rate optimistically degrades exponentially in 
\/L, regardless of the employed purification scheme. It decays, at best, exponentially in L if no 
purification is used. 
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I. INTRODUCTION 



Quantum repeaters enable entanglement distribution between remote parties by relying 
on a network of quantum memory units , [h]]. In their seminal 

paper Briegel et al. demonstrate how to distribute entanglement over arbitrarily long 
distances using ideal quantum memories as well as highly efficient quantum gates for en- 



tanglement purification 



0, 
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131 ] . The resources needed in such a scenario will then only 



grow polynomially with distance for a fixed desired fidelity for the final entangled state. 
Further studies have shown that so long as the coherence time r c of the memories is much 
longer than the transmission delay L/c, where L is the distance between the two parties and 
c is the speed of light in the channel, the above assertion still holds 0, 0] • It is not clear, 
however, how the required resources scale in the limit of large distances when r c is finite. 
Here, we answer this question by assuming that the only error-correction mechanism used 
in the system is entanglement purification in conjunction with allowing for error-free, but 
probabilistic, gates to be used instead of erroneous deterministic ones. We find that, even 



under optimistic assumptions, the system cost explodes as a power of exp [sj (L/ c) 7 T r ] for 
t c <C L/c, unless we employ a fault-tolerant scheme to remedy the memory decay |7J, Il8j ]. 

In order to quantitatively address the cost factor in quantum repeaters, we look at the 
generation rate of maximally entangled states per employed memory in the system. We 
employ proper entanglement measures 14|, instead of merely looking at the fidelity, to 
find this rate. We obtain this rate, in the steady state, assuming that the resources in 
our system are being successively used, according to a proper protocol, to create entangled 
states. Such a rate-over-cost measure is useful for applications in quantum key distribution 
(QKD) ]15], where the generation rate of secure key bits is proportional to the rate of 
entanglement generation. Moreover, it provides us with a fair and practical measure for 
comparing different quantum repeater setups and their contrast with alternative schemes 
for entanglement distribution that do not rely on using quantum memories, such as quantum 
relay structures [16| or the direct transmission of entangled photons. In the latter cases, the 
rate will decay exponentially with distance as a result of loss in the channel. 

Memory decay is one of the most challenging problems in quantum repeater technology 



171 ] . Its deteriorating effect, however, has not yet fully scrutinized. In [7j], authors study the 



role of memory errors in quantum repeaters, but they treat the required initial entanglement 
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as a given resource. This approach cannot fully capture the memory decay problem as it 
neglects to account for the corresponding waiting times during initial entanglement distri- 
bution. Entanglement distribution, regardless of the employed scheme, is a probabilistic 
process, mainly because of its dealing with the loss in the channel, and therefore, the time 
required to entangle two memories is a random variable. Collins et al. 8] consider this 
probabilistic nature in a multiple-memory configuration, and report a numerical-analytical 
rate analysis that includes the effects of memory decay. They model the memory decay by 
associating a lifetime window — within which the stored entanglement is unaffected and be- 
yond which it is destroyed — to each memory. This simple model is not, however, sufficiently 
realistic to properly account for the effect of memory errors on the rate, especially in the 
regime of short coherence times. 

In this paper, we introduce a partial nesting protocol, inspired by the blind repeater 
protocol proposed in [7|], and find the throughput in the two cases of perfect and imperfect 
memories, with or without purification, for a generic quantum repeater system that relies 
on probabilistic schemes for initial entanglement distribution as well as for its entangle- 
ment connection. Such probabilistic architectures for quantum repeaters [2, Q, H] can adapt 
themselves to post-selection-based self-purification mechanisms, which — in the absence of 
memory errors — will allow them to create high-fidelity entangled states over moderately 
long distances without relying on high-quality quantum gates [3j. Such an advantage comes 
at the price of their achieving relatively low entanglement generation rates owing to the 
probabilistic nature of their operation. They are nevertheless attractive options for QKD 
systems whose security over long distances can be guaranteed by a combination of such 
quantum repeater links and sparsely located trusted nodes. To account for the above prac- 
tical issue, in Sec. II, we describe the probabilistic framework of our multi-memory quantum 
repeater system. In Sec. Ill, we present our rate analysis for different scenarios. Numerical 
results will be presented in Sec. IV, and we conclude the paper in Sec. V. 



II. SYSTEM DESCRIPTION 



In our analysis, we study a quantum repeater link of length L with 2™ sublinks of length 
Lq, as shown schematically in Fig. dj with n being the nesting level of the system. At the 
end of each sublink there is a bank of iV quantum memories. We schematically refer to the 
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FIG. 1: (Color online) A quantum repeater with multiple quantum memories per node. At 
each round, we first employ an entanglement distribution protocol to entangle any unentangled 
memory pairs over Dq links. It takes at least Ted = Lq/c to learn about the success/failure of these 
attempts, which succeed with probability P$, hence Ted is the shortest period at which any protocol 
can run. At any such cycle, we also match up entangled pairs at different stations, according to our 
partial nesting protocol (see the text), to perform Bell-state measurements (BSMs), which succeed 
with probability Pm- 

banks in Fig. Q] by nodes 1 to 2 n+1 . Nodes 2% and 2i + 1, for % — 1, . . . , 2" — 1, form physical 
stations Si at which Bell-state measurements (BSMs) take place. For future reference, we 



also denote the family of sublinks of length 2 k L , corresponding to nesting level k — 0, 



, n. 



that connect nodes 2 



k+li 



1 and 2 fc+1 i, for i = 1 2 



n—k 



by D k , and the set of all 
stations corresponding to nesting level k > by S^ = {S 2 k-i(2i~i), 1 = 1,..., 2 n_fc }. 

The probabilistic framework that we use for our rate analysis is as follows. We denote 
the probability of success for our employed entanglement distribution scheme over distance 
Lq by P$. The implicit assumption here is that our entanglement distribution scheme is 
heralding, i.e., there is a mechanism by which we can verify or become informed of whether 
our entangling attempt has failed or succeeded [2|. We also assume that the whole process 
of entanglement distribution from the time that it starts until we learn about its result takes 
Ted = Lq/c. This is the fundamental time period at which we can attempt successively to 
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entangle two specific memories at distance Lq. We further assume that, in the absence of 
memory errors, the memories A and B, upon a successful entangling attempt, are in the 
maximally entangled state \iP + )ab, where |V ;± )a_b = (|01)ab ± |10)ab)/\/2 and \j)k is the 
logical qubit j — 0, 1 for memory = A, 5. Similarly, we associate a success probability 
Pm < 1 to each BSM. This probability, in principle, depends on the state of the memories 
before performing the BSM. In our optimistic analysis, we assume Pm is the maximum 
probability of obtaining a conclusive result for a given BSM module. For instance, if we 
rely on a specific pattern of photodetection events for performing a BSM, the inefficiency 
of photodetectors may result in inconclusive outcomes, in which case we consider the BSM 
has failed. 

The probabilistic nature of our measurement modules imposes stringent conditions on how 
a quantum repeater protocol may run. In the nested purification protocol (NPP) proposed in 
in order to extend the existing entanglement on Dk-i links to links, we must perform 
BSMs at stations S^ k \ for k — 1, . . . , n. Suppose that we have access to deterministic gates 
for local operations and measurements, and that no purification is required. Then, after 
establishing entanglement over all Do links, one can potentially perform all required BSMs 
simultaneously at all stations. In our probabilistic setup, however, there is no guarantee 
that the BSM at one nesting level provides us with entangled states required for the next 
nesting level. More generally, in a multiple-memory configuration, it is important to know 
on which pair of memories a BSM can be performed, and for that one needs to communicate 
between the nodes. This is also the case for the NPP if the employed purification scheme is 
probabilistic 

The above observation requires us to distinguish between two types of BSMs that one can 
perform. In our protocol, we consider a BSM at nesting level k — 1, . . . , n to be "informed" 
if it is performed on a pair of memories in nodes 2i and 2i + 1 known to station Si E S^ k ' to 
be entangled with one of the memories in the relevant distant locations. A BSM is called 
"blind" if the above condition does not necessarily hold 3]. 

With the above considerations, our mth-order partial nesting protocol (m-PNP), for 
m = 1, . . . , n, proceeds as follows. The m-PNP is a cyclic protocol with period Ted- At 
each Tgrj-long cycle, we first attempt to entangle unused memories of D Q links. We learn 
about the success/failure of these attempts by the next cycle. At each cycle, intermediate 
stations also match up as many memories as possible in their two banks for performing 
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(a) 




FIG. 2: (Color online) (a) A quantum repeater link with nesting level n = 2. Here, we have 
only shown memories A2, ■ ■ ■ , A7 that contribute to entanglement generation between memories X 
and Y . In the 2-PNP, all measurements are informed so if measurements labeled by BSM1 occur 
at time t\, BSM2 will occur at time t± +T\. In the 1-PNP, however, all measurements happen 
at time t\. In such a case, BSM2 will be a blind BSM. (b) A snapshot of memories involved in 
creating entanglement between X and Y in the general case of m-PNP. The assumption here is 
that we know that X and ^2 m , ^2 m +i and A 2 m+i, and so on are entangled because all previous 
measurements have been informed. The remaining measurements in 2 n /2 m_1 — 1 stations shown 
above all occur at the same time. Among these BSMs, those measurements that correspond to 
nesting level m, labeled by BSMm are informed, and the rest are blind. 

BSMs. In our protocol, we perform informed BSMs at stations S^ k \ for k — 1, . . . , m, on 
memories known to be entangled over D^-i, and perform blind BSMs on memories known to 
be entangled over D m _i, but not necessarily over Dk-i, at stations S^ k \ for k = m + 1, . . . , n; 
see Fig. [2j These blind BSMs are performed as soon as such measurements can be done. 
Purification may or may not be used over the first m nesting levels. After each measurement, 
or after creating entanglement over distance L, the involved memories will be released to 
be used again in the process of entanglement distribution. The required time for quantum 
measurement /processing is considered to be negligible throughout the paper. 

The key to the polynomial scaling of the rate in quantum repeaters is in performing in- 
formed, rather than blind, BSMs. That can better be understood by comparing the special 

li 

case of m — n, which corresponds to the probabilistic NPP, with quantum relays [l6|. In 
quantum relays, one attempts to establish entanglement over Dq and subsequently performs 
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BSMs at all stations without learning about the success/failure of previous entangling at- 
tempts. Quantum relays can, in principle, be run at a rate faster than 1/Ted and may not 
need to use quantum memories, but their overall rate of entanglement generation will be 
proportional to P^Pm 1 ^ (PsPm) L ^ L ° > which, for PsPm < 1, is exponentially decaying 
with L. This is in contrast with the probabilistic NPP, in which, as we will show, for ideal 
quantum memories, the rate scales polynomially with L as PsPm ^ {L/Lo) log2PM . 

Informed BSMs, however, require us to wait for classical signals. Even if we neglect the 
operation time for local gates/measurements, as we do throughout the paper, it still takes 
Tk = 2 /c ~ 1 Ted, for k = 1, . . . , n, to transfer information from a station in S^ k ' to its closest 
station in S ,< - fc+1 - ) . This implies that there is a minimum delay of ^ED + ^fclTi 1 Tk — T m , in the 
m-PNP, between our initial entangling attempts and the time that we can perform informed 
BSMs at the mth nesting level, during which memory decay is in effect. The memory decay 
will effectively reduce the rate of entanglement generation. An optimum rate may then be 
obtained by a combination of informed BSMs in the first few nesting levels and blind BSMs 
in the remaining final nesting levels. 

To minimize the delay incurred by classical communication and the probabilistic nature 
of entanglement distribution/connection, we have shown that we must use a large number 
of _ per eaeh node of on r qu a„ tum repeater fl Q . To find an optnnistic estate 
of the rate, we then assume that N 3> 1; the rate behavior for a finite number of memories 
has been addressed ,n fi B, Q . 




III. RATE ANALYSIS 

In this section we obtain optimistic estimates for the normalized rate, the rate of entangle- 
ment generation per employed memory, in the two cases of perfect and imperfect memories. 
In the latter case, the only source of error considered is memory dephasing, whose effect on 
the rate is studied in conjunction with different types of purification. 

A. Ideal memories 

Let us first calculate the rate in the case of ideal quantum memories, i.e., when there is 
no source of error, and, therefore, no need for purification in the system. In the m-PNP, 
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starting with no entanglement at time zero, and for N ^> 1, the average number of entangled 
pairs over D n that we create at time T m , is given by = NPsP^ +m ^ ■ Here, NPs 
is the average number of entangled pairs over each D link at time Ted, NPsPJ^" 1 is the 
average number of entangled pairs over each D rn _i link after performing m — 1 levels of 

(2 n — m+1 -m 

informed BSMs, and P^ accounts for the remaining BSMs, where 2 n_m of which 

are informed BSMs for the mth nesting level and the rest are blind; see Fig. 12(b). At time 
Ted, there are on average iV — NPs unused memory pairs in each D link, which can be 
used to create new entangled pairs. We overestimate iV — NPs by N, and repeat the above 
argument to obtain the same value N^ for the average number of entangled pairs created 
over distance L at time Ted + T m . By reusing the same argument for each cycle, in the 
steady state, N^ gives an optimistic estimate of the average number of entangled pairs 
created per cycle. That implies that the steady-state rate of entanglement generation per 
each of N2 n+l ideal memories used in the m-PNP, m — 1, • • • , n, is approximated by 

Ar(") p p (2"- m + 1 +m-2) 

Vm ~ TedA^ 1 2L/c ' 1 ' 

Here, m = n represents the special case of the probabilistic NPP, which is the optimum 
scenario when memories are perfect. 



B. Imperfect memories 

Now, let us consider the case of imperfect memories. Here, we model the memory degra- 
dation by a dephasing process in which a qubit state pa, of memory A, is mapped, after 
decaying for a time period t, to 

r?(p A ) = p(t/2) PA + [1 - p(t/2)}Z A p A Z A , (2) 

where Z A is the Pauli Z-operator acting on memory A, p(t) = [1 + exp(— t/r c )]/2, and r c is 
the memory coherence time. The above dephasing process maps = \i/j ± ) A b{4 ,± \ to the 
following rank-two Bell-diagonal state 

rf ® rf (tt± b ) = P (t)7c± B + [i -v{t)\n B = pi B (t), (3) 

with fidelity p(t). We can also show that a BSM on memories B and C of a four- memory 
system initially in the state p AB (t) ® Pcd(^) l eaves A and D, up to a local unitary, in the 
state p AD (2t). 
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Using the above model, we estimate Rm\ the generation rate of maximally entangled 
states, per employed memory, in the presence of dephasing errors. We first consider the 
no-purification case, and then allow to use purification, without restricting ourselves to any 
specific purification protocol. In each case, for any memories X and Y at distance L in the 
m-PNP, we obtain the final state p^xy upon establishment of entanglement between them. 
We denote by Aj the memory in node i = 2, . . . , 2 n+1 — 1 that contributes to the successful 
creation of entanglement between X and Y; see Fig. [2j 

To underestimate the decay effect, we make a further assumption in calculating p/ t xY- 
In the m-PNP, after performing the last set of BSMs, corresponding to nesting levels m and 
higher, we must in principle wait, before claiming that X and Y are entangled, until we 
receive the classical data from all middle nodes. The effective state in such a case must then 
include the decay during this last step of communication. In certain applications such as 
QKD, however, we can perform the required QKD measurements on X and Y concurrent 
with the time that blind BSMs are performed, and later confirm whether they are entangled 
or not. Under the latter assumption, we do not need to account for the classical delays that 
occur after the last set of BSMs corresponding to nesting levels m and higher. Throughout 
the rest of this Section, we assume that we are operating under such conditions. 

With the above assumptions, we overestimate the number of maximally entangled states 
that can be obtained out of Nm^ copies of Pf,xY by Nm^ E c (pf )X Y)- For a state p, E c (p) 
is the entanglement cost, viz. the ratio M/K for starting with ^+)« ,M an d obtaining 
p® K , in the limit of large K, using local operations and classical communication (LOCC). 
It is known that for the asymptotic version E°°(p) of any "well-behaved" entanglement 
measure, we have -Ed(p) < E°°(p) < Ec(p), where -E'd(p) is the distillable entanglement, 
i.e., the ratio M/K for distilling jP K to (/fr+)® M \ in the limit of large K, using LOCC 



14| . Hence, our choice of entanglement cost as our entanglement measure is in line with 



our other optimistic assumptions. For pjy(t), we have E D {p^ Y ii)) — 1 — H{p{t)) and 



Ec(p XY (t)) = H{\ + y/p(t)[l-p{t)]), where H(p) = -p\og %P - (1 - p) log 2 (l - p) [l4|. 
Here, Ep(t) also gives the asymptotic yield for the one-way hashing protocol for purifying 

P+y(t) |l2]- 
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1. Without Purification 



Let us first consider the case where no purification is employed. Suppose that we attempt 
to entangle D links at time zero. Then by the time that we become informed of the estab- 
lishment of entanglement, memories have already decayed for a period Ted. For instance, X 
and A2 are in the state Px^ 2 (Ted)- At Ted = Ti, we also perform the BSMs corresponding 
to the first nesting level. That leaves X and v4 4 , right after performing the BSM and up to 
a local unitary, in /5j^ 4 (2Ti). It takes T\ to inform A4 of the BSM success, hence at time 
T2 = T E d +Ti, X and A 4 are in p^ A (3Ti). By continuing the same argument, we can show 
that at time T m , p/,xy = PxY^'m), where t^m = T n+1 + (m— l)T n , for m = 1, . . . , n. By time 
T m , we do not necessarily know whether X and Y are entangled, but as mentioned before, 
in some applications such as QKD, we can often use X and Y at time T m , and later verify 
whether they are entangled or not. Under this assumption, our optimistic rate is given by 



Given that, to the first order, Ec{p^ Y ^)) °^ £exp(— 2t/r c ) for t 3> r c , and that tm > L/c, 
the rate in Eq. (jl]) exponentially decays with distance for a finite value of r c . This implies that 
the partial nesting protocol cannot help us avoid the exponential decay of the normalized 
rate with distance if we do not use any purification. 

2. With Purification 

Now, let us consider the case when we allow to use purification. In the m-PNP, our last 
chance to purify memories is in the mth nesting level. In order to obtain an upper bound 
on the rate, suppose that at the nesting level m — 1 > 1, we are provided with maximally 
entangled states tx + over -D m _ 2 links. (We consider the special case of m — 1 later.) The 
BSMs at stations 5 , ^ m_1 \ upon success, will leave us, up to a local unitary, in tt + . There 
will be, however, a delay of T m -\ to inform stations of the success/failure of the above 
BSMs, during which the memories decay to p + (T m _i). In the special case of m = 1, the 
minimum delay is given by Ted, hence the decayed state will be /3 + (Ted)- The rate is then 
optimistically approximated by 



Bff « Q^E c (p+ Y (W)), (no purification). 



(4) 



« Q$E c (p+ Y (TW)), (with purification), 



(5) 
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where = max(TED, T m -i). Here, for a fixed value of n, Ec(p\ Y (Tm^)) is a decreasing 
function of m, whereas Q$ is an increasing function of m. That implies that, for fixed 
values of r c and L, there are optimum values of m and n that maximize the normalized rate, 
and such a maximum may scale better than exponentially with distance, as we show next. 

In order to find the asymptotic behavior of R opt = max m n R$ in the limit of large 
distances, we consider all possible choices that we have for our system parameters m and n. 
The only additional assumption in our analysis is that P$ oc exp(— aL ), where a is inversely 
proportional to the attenuation length of the channel. This is a reasonable assumption as 
most proposed entanglement distribution schemes require the transmission of a single photon 
along the channel 2, 5, f|. Now let us consider the following cases for the optimum values 
of m and n, respectively, denoted by m op t and n op t, all in the limit of large distances: 

(1) Suppose n opt is finite, then, because Pg oc exp(— aL/2 n ), R opt degrades exponentially 
with L. 

(2) Suppose n opt — > oo but n opt — m opt is finite. Then the entanglement cost term in Eq. (j3j), 
will scale as exp[— L/(cr c 2 n ° pt ^ mopt+1 )], which again represents an exponential decay of the 
normalized rate with distance. 

(3) Suppose n opt — > oo and n opt — m opt — > oo, then 

R opt oc exp(-«Lr ) exp ( ^^ j E c (p+ Y (2— 2 L°^ /(cr c ))), (6) 

where L° pt = L/2 n ° pt . (Note that R^ < R { ?\ therefore T^ = 2 m °^- 2 If^ / 'c.) Now, if 
2 m ° pt LQ Pt is finite, then the entanglement cost term in the above equation will be a constant, 
whereas the second term, so long as Pm < 1, decays exponentially with L. If 2 m ° pt L°^' t — > oo, 
we can then replace the entanglement cost term in Eq. ([6]) with its asymptotic value to obtain 

R opt oc exp(-«Lr ) exp ( ^f^t ) exp(-2— 1 (7) 

In the above equation, there are two competing terms: the second term that decays exponen- 
tially with L/(2 m °p t_1 Lo Pt ) and the third term that decays exponentially with 2 m °p t ~ 1 Lo Pt . 
The best rate-versus-distance scaling can be achieved when these two terms scale the same 
with L. That can be achieved by assuming 

2 m op t-i iv opt = for 

some (5 > 0, in which 

case both above-mentioned terms decay exponentially with a/L. After optimizing over (3 
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and assuming that m op t — > oo, we find that 



. L\n(l/P M ) \ 
R opt oc exp | -2 J , (8) 

which is exponentially decaying with y/L, rather than L as obtained for cases (1) and (2). 
This is in fact the best rate-versus-distance scaling achievable for our quantum repeater 
system under our optimistic assumptions. Given that the above optimized rate is obtained 
in the regime of m p t ,n opt — > oo, we can further manipulate the ratios R^+i I Rm and 
Rm+i/R*n to find that, for large distances, the optimum rate is obtained at 

«21nP^/«, (9) 
which is a constant identical to what one obtains for the ideal-memory case, and 



2 m ° pt « a\j2LcT c / In P^ 1 . (10) 

This shows that the optimum rate is obtained when 2™ opt oc L and 2 m ° pt oc yl- Such 
an improvement over the no-purification case is associated with the employed purification 
scheme as we discuss next. 

The proposed purification schemes can be divided into one-way and two-way protocols 



12J. In most proposed one-way protocols, it is commonly the case that after performing 
proper local operations on one side, the measurement results are sent to the other side at 
which, after further manipulation, the purified states are obtained. In other words, the 
operation of such schemes is equivalent to that of a quantum code, in which errors occurred 
during both the processing as well as the storage time may be corrected. In most proposed 
two-way purification protocols, it is commonly the case that we perform local operations on 
both sides, and exchange the measurement results between the two parties. The purification 
process can then be completed only after the receipt of classical data. The problem with 
these protocols is that the faulty memories will decay during the transmission of classical 
data. That may result in an exponential decay of rate with distance as detailed below. 

To see the effect of a typical two-way purification scheme on the normalized rate, suppose 
that in order to purify pab, an arbitrary two-qubit state between two partners A and B at 
distance I, we apply local operations/measurements on M copies of pab to obtain K < M 
copies of p A U g. Then, we communicate between A and B, to verify whether the purification is 
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successful. Upon success we perform another set of local operations, modeled by a quantum 
operation A, on (T^ c <g> Yf/J® K [(p A u B )® K }. For many schemes proposed in [ll], A is just the 
identity operator due to a discard-or-keep action, so we further assume that the dephasing 
process modeled by ® Yf, commutes with A. Under these conditions, we show that the 
fidelity of memories A and B after purification cannot exceed p(l/c). Suppose p AB = A(p A u B ), 
where F = tr(it + p AB ) is the fidelity of the state after purification without considering the 
memory decay The fidelity of such state after memory decay for a period t = l/c will then 
satisfy 

tr^rf ® rf {p F AB )\ = tr[rf ®Yf{^ AB )f AB \ 

= tT[(p(t)nX B + [l-p(t)]n AB )p F B } 

< p{t)F + [l-p(t)){l-F)<p(t), (11) 

where the first equality comes from the cycling property of the trace operator, and equality 
holds for the first inequality if p AB = Ffr AB + (1 — F)n AB . We used F < 1 in the last 
inequality This implies that even if p v AB = tt ab , by the end of the above purification 
procedure, we ideally end up with paB + (1/c). Using such purification schemes in the m- 
PNP will then leave memories X and A2m, at best, in the state P^ j4am (T m ), right before 
performing the last set of BSMs. This is the case too for other similar pairs of memories. 
After performing the remaining BSMs, we end up with Pf t xY — Pxy(^n+i)> f° r which the 
entanglement cost scales as exp[— 2L/(ct c )], and that proves our claim. 

Whereas the two-way purification schemes considered above cannot provide us with 
the optimum rate-versus-distance scaling shown in Eq. ([8]), there are one-way purification 
schemes that achieve a similar scaling. For instance, in the m-PNP, suppose we perform no 
purification up until the m-th nesting level, at which point we perform a one-way hashing 
protocol [12I on entangled states distributed over -D. m _i. These entangled memories, before 
applying the purification and up to a local unitary, are in the state p + (tffl /2). To perform 
the hashing protocol, we need to perform encoding measurements on one side of each entan- 
gled link over -D m _i and send the results to the nodes on the other side of the link. Suppose 
the nodes that receive this encoding information are those at which the mth nesting level 
measurements will occur. For all other nodes, we not only can perform, at time T m , the mea- 
surements required for the hashing protocol but also, at the same time, we can perform those 
blind BSMs, according to our partial nesting protocol, needed for entanglement swapping as 
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well as the required measurements on the end nodes with regard to our QKD application. By 
this trick, only memories located at will undergo an additional T m -long decay due to the 
transmission delay before getting purified. The normalized rate after performing the above 
purification as well as the remaining BSMs is then given by Qm E D {p^ Y (tm + T m /2)), 
which, for m = 1 and L oc \/L, degrades exponentially with \fL similar to Eq. (jBJ). Here, 
we used the fact the asymptotic yield of the one-way hashing protocol for our decayed state 
is given by the entanglement of distillation E D (p~x Y (t)) oc exp(— 2t/r c ), for t 3> r c . 

The above example shows that with even one step of purification we can achieve the 
optimum scaling achievable by a quantum repeater system that only relies on entanglement 
purification for mitigating the effect of memory errors. This is a weaker requirement than 
what needed for fully fault-tolerant schemes, which rely on the frequent use of error correction 
schemes. Although here we assumed that the gates used in the purification scheme are 
deterministic and error-free, we can possibly relax this assumption by using probabilistic 
gates for purification. How the rate decays in this new scenario is a matter of further 
investigation. 



IV. NUMERICAL RESULTS 



In this section, we numerically study the effect of different system parameters on the 
normalized rate. In all graphs presented, we assume Ps = 0.2 x io~ aoiLo ; which is an 
appropriate choice for the schemes proposed in . Furthermore, The underlying physical 
channel is assumed to be an optical fiber with 0.17dB/km loss, for which c = 2x 10 8 m/s and 
1/a = 50 km. For numerical purposes, we have used the normalized rate given by Eqs. (j4]) 
and (jSJ), for, respectively, systems without and with purification. 

Figure EJ^a) shows the rate in Eq. (jHJ) versus nesting level n for L = 1000 km, hence 
L/c = 5 ms. For each value of n, we have optimized m to maximize the rate. The optimum 
values of m turned out to be equal to n for all the marked points in Fig. [3](a). This figure 
clearly shows the criticality of L/c as a benchmark for the required coherence time. Whereas 
the rate drops over an order of magnitude when the coherence time is reduced from 5 ms 
to 1 ms, there is only a slight improvement in the rate if we use memories with 100 ms 
coherence time. We should bear in mind, however, that these values are only valid in the 
case of N ^> 1. The benefit of using quantum memories with longer coherence time is that 



14 




3 4 
nesting level, n 



0.01ms 0.1ms 1ms 10ms 100ms 
coherence time, r c 



Is 



10s 



FIG. 3: (Color online) An optimistic estimate to the generation rate of maximally entangled 

(n) 

states per employed memory, R m , for the m-th order partial nesting protocol, in the presence 
of memory dephasing. The implicit assumption is that we employ a sufficiently large number of 
memories to minimize the waiting times due to classical communication. In (a), we have plotted 
B&\ optimized over m, versus the nesting level, for three different values of coherence time. In (b), 

(n) 

we have plotted Rm versus coherence time, optimized over m and n. For fixed values of r c , the rate 
scales exponentially in y/L. The pairs (n, m), on several points on the graph, show the optimum 
nesting levels, n, and the optimum numbers of informed BSMs, m, at those points. In both (a) and 
(b), we assume that our entanglement distribution succeeds with probability Ps = 0.2 x io~°- 01L ° 
and the total distance is L = 1000 km, thus L/c = 5 ms in optical fibers. 



we can still achieve the same value of P^m with a lower number of employed memories 10 ] 
Another interesting point in Fig. E^a) is that the achievable rate for both values of n = 3 
and n = 4 is about the same. Although, numerically, n = 4 attains the maximum rate, for 
practical purposes, it is beneficial to use the less costly architecture corresponding to n = 3. 
One last point about this figure is the difference in the slope of rate for values of n smaller 
than the optimum value of n, and those that are larger. It is because for n < n opt , loss is the 
dominant factor, whereas for n > n opt the BSM success probability brings the rate down. 

Figure 0(b) shows the optimum rate versus coherence time for three values of Pm at 
L = 1000 km. On several points on the graph, we have indicated pairs (n, m) that achieve 
the maximum achievable rate. For r c L/c, the partial nesting protocol favors using blind 
BSMs as they require memories with shorter coherence times. For lower values of Pm, one 
needs to use fewer nesting levels to attain the maximum achievable rate. The optimum 
values of n in the above cases agree with the relation Lg pt w 2\n(P M 1 )/a found in the 
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FIG. 4: (Color online) Estimates for the generation rate of maximally entangled states per 

(n) 

employed memory, Rm j versus coherence time, optimized over tn and n, for Pm — 0.75, 
Ps = 0.2 x io~°- oli o, = 1000 km, and L/c = 5 ms. The employed entanglement measures 
in these graphs are entanglement cost for the solid lines, and entanglement of distillation for the 
dashed lines. For fixed values of r c , the curves scale exponentially in L if no purification is used and 
exponentially in \f~L with purification; see dash-dotted line for comparison. The two-tuples (n, m), 
on several points on the graph, show the optimum nesting levels, n, and the optimum numbers of 
informed BSMs, m, at those points. 



previous section. 

Figure H] demonstrates the rate behavior in both limits of long and short coherence times. 
The solid lines represent the optimistic normalized rate when entanglement cost is used as 
our entanglement measure as in Eqs. (j3J) and (jSJ). The dashed lines represent the achievable 
normalized rate when entanglement of distillation is used instead. It can be seen that the 
difference between the two cases is small. It is noticeable that the normalized rate, regardless 
of the employed entanglement measure or purification scheme, approaches a constant value, 



maXnQ^, when r c ^> L/c, in accord with what reported 



in 0, 



and Eqs. (ED and (jSJ). In 
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the other unexplored extreme of L 3> cr c , however, the situation is much different. In the 
no-purification case, because, to the first order, Ec(pxyi^)) ^ texp(—2t/r c ) for t ^> r c , the 
rate in Eq. (j3J exponentially decays with distance. Proper use of purification helps, but only 
a little. As shown by the dash-dotted line in Fig. HI at the optimum values of m and n, 
in Eq. §5§ scales as a power of exp[— a/ (L/c)/t c ]. The same effect has been shown in Fig. El 
where we have plotted the optimized rate versus L for several values of coherence time. It 
can be seen that for r c = 100 ms, the logarithm of rate scales linearly with logL, which 
represent a polynomial scaling of rate with distance. For short coherence times, however, we 
will observe the exponential decay with a/Z as discussed before. Another interesting point 
is the optimum values of m and n, which scale, respectively with log a/Z and log L. 



V. CONCLUSIONS 



In this paper, we studied quantum repeaters that relied on probabilistic schemes for 
entanglement distribution and connection. We introduced a proper rate-over-cost figure of 
merit for such systems by looking at the ratio between the rate of entanglement generation — 
when the system resources are being fully employed to successively create entangled states— 
and the total number of employed memories. We used entanglement measures to quantify 
the generation rate of entanglement. We believe that to properly compare different quantum 
repeater setups, one must study such a normalized rate instead of merely considering the 
fidelity of generated entangled states at the end of only one round of the protocol. We studied 
the effect of memory dephasing under the assumption that the only mechanism employed to 
mitigate errors is purification. We further assumed that the number of memories employed 
per node is sufficiently large such that, by parallel entangling attempts, we can minimize the 
unavoidable communication delay between relevant nodes. The memory decay, in principle, 
can be dealt with using complex fault-tolerant schemes [l8], which may not be available 
in the near future, and has not been considered in our work. We showed that, whereas 
in the ideal case of infinitely long coherence times the normalized rate scales polynomially 
with distance, with imperfect quantum memories the rate-over-cost asymptotically scales, 
at best, as a power of exp [— a/(L/c) /r c ] , with r c being the memories' coherence time. This 
rate behavior is nevertheless superior to the fully exponential decay of rate with distance 
in quantum relays, which do not use quantum memories, or in systems that do not employ 
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FIG. 5: (Color online) Estimates for the generation rate of maximally entangled states per em- 

(n) 

ployed memory, Rm , versus total distance L, optimized over m and n, for different values of 
coherence time at Pm = 0.5 and P$ = 0.2 x io~°- oli °. The employed entanglement measure in 
these graphs is the entanglement cost and a proper purification protocol is used. It can be seen that, 
for t c 3> L/c, the curves scale polynomially with L, whereas, for r c -C L/c, they scale exponentially 
with yL. The two-tuples (n, m), on several points on the graph, show the optimum nesting levels, 
n, and the optimum numbers of informed BSMs, m, at those points. They, respectively, scale with 
log L and log \/L, for large distances. 

purification but use decaying memories. Our analysis helps us estimate the minimal cost 
of our quantum communication systems. In the particular case of QKD, to create 1000 
secure key bits/s over 1000 km, we may need thousands of memories with coherence times 
exceeding or on the order of 10 ms. 
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